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Abstract: We construct a Q = 1 supersymmetry and U{1)^ global symmetry preserving 
deformation of the type IIB matrix model. This model, without orbifold projection, serves as a 
nonperturbative regularization for Af = 4 supersymmetric Yang-Mills theory in four Euclidean 
dimensions. Upon deformation, the eigenvalues of the bosonic matrices are forced to reside 
on the surface of a hypertorus. We explicitly show the relation between the noncommutative 
moduli space of the deformed matrix theory and the Brillouin zone of the emergent lattice 
theory. This observation makes the transmutation of the moduli space into the base space 
of target field theory clearer. The lattice theory is slightly nonlocal, however the nonlocality 
is suppressed by the lattice spacing. In the classical continuum limit, we recover the = 4 
SYM theory. We also discuss the result in terms of D-branes and interpret it as collective 
excitations of D(-l) branes forming D3 branes. 
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1. Introduction 

The study of supersymmetric deformations of the type IIB matrix model (or IKKT matrix 
model [1]) is interesting both from the field theory and string theory points of view. The de- 
formations may alter the classical moduli spaces into noncommutative spaces, and hence, may 
form a bridge between zero dimensional, finite rank matrix theories and higher dimensional 
regularized quantum field theories. In this work, we show such a correspondence between a 
Q = 1 supersymmetry and U{1)^ global symmetry preserving deformation of type IIB matrix 
model, and = 4 supersymmetric Yang-Mills theory in four dimensions. Mainly the de- 
formed type IIB matrix model serves as finite Nc matrix regularization (which will be shown 
to be equivalent to a lattice regularization), where Nc is the number of colors, for = 4 
supersymmetric Yang-Mills theory in four dimensions. This construction can be thought as 
supersymmetric generalization of the relation between twisted Eguchi-Kawai (TEK) and four 
dimensional lattice gauge theory [2,3]. 
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The essential point of the construction is the relation between two seemingly different 
concepts. The noncommutative moduli space of vacua, which is the zero action configuration 
of the deformed bosonic action, and the noncommutative base space, the space that the 
target field theory lives in. In the type IIB theory without deformation, the eigenvalues of 
the bosonic matrices can be anywhere in the moduli space. The deformations, that will be 
considered in this paper, alter the structure of moduli space in such a way that the eigenvalues 
of bosonic matrices are forced to reside on the surface of a four dimensional hypertorus, T^. 
This is essential for the emergence of four dimensional spacetime and will be discussed in 
detail in sections |^ and 4.1. 



In the description of the deformed type IIB action in section |2|, we use the notation 
introduced in [4], which makes the SU{5) x U{1) subgroup of global SO{W) R-symmetry 
and Q = 1 supersymmetry manifest. Therefore, we use five complex matrices and their 
hermitian conjugates Zm as the bosonic degrees of freedom instead of ten bosonic hermitian 
matrices. This notation is more akin to our four dimensional habits and the bosonic action 
becomes just a sum of D-term and F-terms, which are derivatives of superpotential. The 
bosonic action of the deformed type IIB matrix model is given by 



osonic — 9 

9^ 



J*,nn/2 - - 1 2 



a.i) 



where the antisymmetric phase matrix ^mn are the deformation parameters. This change 
in the action may be considered as the deformation of the "superpotential". The action 
Eq. (1.1) may be obtained by the dimensional reduction of the four dimensional hypercubic 
lattice action of = 4 SYM to zero volume (a point) upon the use of twisted boundary 
conditions. The phases ^^u, where /i, = 1, . . . , 4, are associated with electric and magnetic 
background fluxes [5] and <I>^5 = — X]a=i*^Ma- We will show that a two parameter family 
of deformation which, for example, corresponds to electric and magnetic fields in z-direction, 
suffices to generate the target M = 4 theory ^ . 

The technique described here is different from the recent constructions of lattice regular- 
ization for extended super symmetric Yang-Mills theories [4,11-13] and deconstruction of six 
dimensional theories [14,15]. The current approach does not require orbifold projection of a 
mother theory. As it will be shown explicitly, the deformation of the mother theory becomes 
the lattice regularization. Recall that the orbifold projected matrices are sparse and encode a 



^This type of deformation appeared in the EK context under the rubric "twisting" [2,6,7], and in the 
context of supersymmetric gauge theories in four dimensions as a marginal deformation of superpotential (as 
l3 deformation) of A/" = 4 theory to A/" = 1 (see for example [8-10]) 
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commutative spacetime lattice. In our approach the spacetime lattice is generated by the zero 
action configurations of the deformed type IIB matrix model and there are no constraints on 
our matrix variables. Similarly, based on [11] and [3,16], the lattice action for noncommu- 
tative counterparts of the extended supersymmetric gauge theories are also constructed, by 
using orbifolds with discrete torsion, in [17] for spatial and [18] for Euclidean two dimensional 
lattices. (For recent progress and independent approaches on supersymmetric gauge theories 
on lattice, also see [19,20]. For a review on discrete torsion, see [21] and the references in [18]) 

To show the equivalence of the deformed type IIB matrix model to the four dimensional 
theory on a discrete torus, we decompose the bosonic matrices into a background and fluc- 
tuations around it. By using the background field techniques, it is possible to exhibit that 
the fluctuations around the zero action configurations can be identified as a lattice action 
on a four dimensional discrete torus. This action has many parallels to the lattice action in 
ref. [4]. The main difference with respect to the action in [4] is the noncommutativity of the 
base space lattice, which can be traded with the nonlocalness of the interactions. However, 
we show that two types of finite-volume continuum limits are possible. In a particular scal- 
ing of the deformation parameter (in which the noncommutativity length tends to zero with 
respect to the linear extent of the box), the interaction localizes in the position space and 
the outcome is ordinary local field theory. In the limit where the noncommutativity length is 
nonvanishing, the outcome is nonlocal AA = 4 SYM. 

It is also possible to interpret the result of this paper within string theory. The main result 
can be understood as a collective excitation of a finite number of D(-l)-branes to generate 
discretized D3-brane(s) wrapped on a four dimensional hypertorus. Notice that we do not 
use compactification (see [22] for a review) to generate the higher dimensional theory. Also, 
our D(-l)-brane configurations are inherently different from the BPS-saturated backgrounds 
discussed in [23] , which have typically finite action and are unstable at finite N^- We discuss 
this direction briefly in section |^. The rest of the paper is organized as follows: In section ^ 
and appendix |A|, we set the notation and derive the deformed type IIB action. In section ^, 
we analyze the noncommutative moduli space and make the relation with TEK more explicit. 
In section ^, we show that the the deformed type IIB action can be rewritten as a lattice 
action and obtain the classical continuum limit. 

2. The type IIB matrix model and its supersymmetric deformations 

The type IIB matrix model can be obtained by the dimensional reduction of the ten dimen- 
sional M = \ U {Nc) supersymmetric Yang-Mills theory to zero dimension. The theory possess 
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a Gr = 50(10) global R-symmetry, inherited by the (Euclidean) Lorentz symmetry of the 
ten dimensional theory. The bosonic and fermionic degrees of freedom of the matrix theory 
are in vector (10) and positive chirality spinor (16) representations of 5*0(10). The action 
of the type IIB matrix model is conventionally expressed as 

10 

5 = -^ V -Ti [v^,vpf+'-TTuo'^CT^[v^,uo\ (2.1) 

where Va = v^t^ is a Hermitian matrix and lo = uj"'t^ is a 32-component spinor which 
has only its upper sixteen component nonvanishing. The ten gamma matrices form a 
chiral basis of 50(10) and C is charge conjugation matrix. The Lie algebra generators are 
normalized such that Tr f^t^ = 5°'^ . This action respects Q = 16 supersymmetries generated 
by 5vci = —K^CTaU), and 5uj = j [wa, f/j] [Tq,, r^]fi; , where k is a positive chirality Grassmann 
spinor. 

We aim to construct deformations of the type IIB matrix model, which preserves Q = 1 
supersymmetry and U{1)^ global symmetry. Since the bosonic variables that appear in the 
action Eq. ( |2.1D are not eigenstates of the C/(l)^, this form of the action is not adequate for 
our purpose. The type IIB action, in terms of eigenstates of SU{5) x U{1), which provides 
the right basis for our purpose, is constructed in [4] and is outlined in the next section. 

The motivation behind this construction is to obtain a finite rank matrix regularization 
of = 4 SYM theory in four dimensions. It is well-known that the TEK model is a matrix 
regularization of (noncommutative) Yang-Mills theories in even dimensions. Our goal here is 
to construct a fully supersymmetric counterpart of this relation. It is, therefore, important 
to construct the supersymmetric matrix model analog of TEK. The TEK is obtained by 
dimensional reduction of a fully regulated lattice gauge theory to a single site by using the 
twisted boundary conditions. We will follow the same strategy and reduce the lattice regulated 
= 4 SYM to a single point. 

2.1 M = A SYM on hypercubic lattice 

In order to construct the deformed Type IIB matrix model, we start with the four dimensional 
hypercubic orbifold lattice. In [4], the lattice action of = 4 SYM is given for both A'^ 
lattice and hypercubic lattice. Here, we will use hypercubic for convenience. The dimensional 
reduction of the lattice to a single site by using periodic boundary conditions reproduces the 
type IIB matrix model in the eigenbasis of U{1)^ symmetry. The reduction by using twisted 
boundary conditions gives the deformed Type IIB action. 

Let us first rewrite the action of the = 4 SYM theory in d = 4 dimensional hypercubic 
lattice and establish our notation. The fundamental cell of a hypercubic lattice contains one 
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fields 
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0-cell 
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l-cell 






2-cell 






3-cell 




-\- Gp -\- Gfj Gp G5 


4-cell 




65 — 2^^=i 65 



Table 1: The distribution of the bosonic and fermionic helds to the p-cells of the hypercubic lattice. 
Each field is totally antisymmetric in its indices. For convenience, we will also use ^'^'"^ — e'^''°'^^^5, 
{Zfj.upa, i'fj.upa) — £pvpa{z^, V'^): s,nd z^^''f"^ = £'^''^'^25. The cartesian basis vectors are given by (e^)^ = 
(5^1/ where fi^v — 1, . . .A . 

site, four links, six faces, four cubes and one hypercube, collectively called p-cells. There are 
bosonic and fermionic fields associated with each p-cell. For a d = 4 dimensional lattice, the 
total number of p-cells in a unit-cell is 

4 

2^ = ^Qj = le4e6©4©l (2.2) 

The fermions, which form the 16 dimensional spinor representation of the 50(10) are dis- 
tributed among the p-cells. The labeling of the fermions follows the name of cell that it lives 
in: w ~ 16 = A © © ipu © i^""^ © ippvpo- = 1©4©6©4©1 where each fermion is totally 
antisymmetric tensor and they reside on sites, links, diagonal of faces, diagonal of the cubes, 
and the body diagonal of the hypercube [24] The upper and lower indices determine the ori- 
entation of the fermionic link field (see Table.|^). For notational convenience and brevity, we 
introduce = ^-H^^^P'^ and tp^ = It^ipf^'^P'^ 

The bosons, which form a 10 dimensional vector representation of SO{10) are associated 
with 1-cells and body diagonal of 4-cells. We have ~ 10 = © © z^i/po- © z^^^^" = 
4 © 4 © 1 © 1. For convenience, we rename z^ = ^^^^z^^pa- and Z5 = Ie^z'^'^p^ 3_ 

Naturally, the Q = 1 superfields on the hypercubic lattice are associated with p-cells. 

^The idea of associating the p-cells with p-form fermions also appears as a natural implementation of 
Dirac-Kahler fermions on the lattice. See for example [25] and reference therein. 

^Throughout this paper, we adopt the convention that repeated indices are summed, where a,(3... are 
SO(IO) indices summed over 1, . . . , 10; the indices pL,v, . . . are SU(4) indices summed over 1, . . . , 4; m,n, . . . 
are SU (5) indices summed over 1, . . . , 5. In the description of the continuum (which will be used only in section 
4.2 and 4.3 ) the indices a,b, . . . are 50(6) _R-symmetry indices summed over 1, . . . , 6; and . . are S'0(4) 

spacetime indices summed over 1, ... ,4. We hope the latter wont cause an inconvenience for the reader. 
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The component form of the superfields and the p-cell that they are residing in are given as 



0- 


cell : 


A„ = A„ - iOdn , 


1 - 


cells : 




2- 


cells : 




3- 


cells : 


S/i5,n = ^/i5,n — 29 E 


4- 


cell : 


7l = zi + ^/29i,i , 



(2.3) 



where the z^^^i and 25 „ are supersymmetry singlets and the dn is an auxiliary field. The 
antiholomorphic functions E^,^^^ and E^^^^ are given by 



E 



(2.4) 



The antiholomorphic E functions are supersymmetry singlets as well since they only depend 
on z, but not z. They can be thought as zero dimensional counterpart of the superpotential 
(more precisely, the F-term) 

The hypercubic lattice action may be written in manifestly Q = 1 supersymmetric form 
as in [4] ^ 

^ ~ "^2" ^ ^ J ^^'^■^ndo-^n A.n(^m,n— em^m,n— em ^m,iiZm,n) + 2 

V2 " 



+ 



mnpqr 



^mn,n{Zp,n-ep^qr,n+ein+en ^qr,n-eq-er^p,n+em+en) 



(2.5) 



where the holomorphic E functions are 



En 



Zn,nZn^,n+en) • 



(2.6) 



The E function is a function of superfields Z and has a 9 expansion where the lowest compo- 
nent of E is the minus hermitian conjugate of E: 



E = £; + edraEi}"', E 



(2.7) 



"* The description of plaquette types and the flux passing through them are easier to understand when we 
make the SU{4) indices iJb,v . . . manifest, However, it is more economical to use the SU{h) notation occasionaUy 
(by keeping in mind the geometric structures that the fields are associated with). Hence, wc will continue to 
use both. Hereafter, we sometimes combine the ^lupa index which runs over 1, ... 4 with 5 and use the Latin 
letters mnpqr which ranges over 1, ... 5. Mainly, we combine the four link multiplets and one 4-cell multiplet 
as + -> Z"" and 5^ + 25 — > Zm, the six fermi multiplets on the faces and four fermi multiplets on the 
3-cells combine into H,^i. + H^s — > Smn and naturally iS^jy + S^s Emn- 
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Notice that the last term in the action is not integrated over the superspace, yet it is respectful 
to Q = 1 supersymmetry. It is not hard to show that the 6 component of the last term is 
identically zero due to Jacobi identity, and translational invariance of the lattice action. Hence 
it is supersymmetric. 

The dimensional reduction of the lattice action Eq. ( p.5| ) to a single site by imposing the 
periodic boundary conditions (for example A^+g^ = An) and by integrating out the auxiliary 
d field yields the type IIB matrix theory action Eq. (|2.l| ) expressed in terms of SU (5) x U (1) 
multiplets 



9' 



^ (2. 



'2 

m=l m,n=l 



1 \/2 



This action is merely a rewriting of Eq. (2J), in terms of a new set of variables. Recall that 



Eq. (|2.lD is expressed in terms of SO{10) multiplets, vector 10 and positive chirality spinor 



16. We could have obtained Eq. (|2.8| ) from Eq. by just employing the decompositions 
~ 10 ^ z"" e ~ 5 © 5 for bosons, and ~ 16 ^ A © -0™ © ~ 1 © 5 © 10 for 
fermions into SU{5) multiplets as in ref. [4]. The dimensional reduction is used to set the 
notation and for later convenience. 

Notice that the bosonic matrix model takes a form similar to the potential in d = 4 
dimensional M = 1 supersymmetric theories. There the potential is a sum of d and F 
(derivative of the superpotential) term contributions in the ^ + \F\'^ form. In the zero 
dimensional matrix model expressed in manifest Q = 1 superfields, the role of F is replaced 
with the E functions and the bosonic action can be written as 

^2 _ ^ ^ - 2 

Sbosonic l~ ^d [Zj^i, Z ] -|- ^ ^ I-S'mni (2'9) 

where E^™" = [ z"^ , z" ] . The deformation that we will construct in the next section will alter 
E, but not the d term. 

2.2 Dimensional reduction and twisted boundary conditions 

To formulate the deformed type IIB model, we start with the lattice action Eq. (|2.5| ) of the 
= 4 SYM theory in four dimensional hypercube. Then we reduce the lattice theory to a 
single site by imposing the twisted boundary conditions. This is indeed the same philosophy 
used in the construction of the twisted Eguchi-Kawai model [2] . 
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The difference between deformed Type IIB and TEK is that in = 4 lattice action, there 
is a superfield associated with each p-cell. Hence, there are many gauge invariant plaquettes 
terms possible, which are shown in figj2| and fig.|l], and will be explained momentarily. In 
TEK, we only have unitary link (1-cell) fields and the only possible gauge invariant plaquette 
terms are the square plaquette terms. 

In reducing the lattice action Eq. ( |2.5D to a single point by using twisted boundary 
conditions, we request periodicity modulo gauge rotations g^. Let ^' be a a generic field 
associated with a p-cell on the lattice. We utilize 

4 

^n+es = 55^n55' where g^ = {Y{gu)\ (2.10) 

1^=1 

and then set to a site independent value. The global gauge rotation matrices 3^ satisfy 
the 't Hooft algebra [5] 

g^^g,. = e"^^''gugi„ equivalently, e'^'^^^^/'^g^g^ = e'*'"'/^^^^^ (2.II) 

The antisymmetric matrix is identified as the flux passing through the fiv plane. More 
precisely, = 2tt^-^ where n^p is an integer associated with the integer electric (n4j) 
or magnetic (riij) flux modulo A'c. Also, note that ^'n-e„ = gt^ngu- The second form of 



the equation Eq. (2.11) is symmetric under the exchange fx ^ u and will be more useful in 



appendix. The commutation relation of g'^ with g^ is given by 

4 

5/^55 = e'*''^559/., where $^5 = - ^ (2.12) 

The $^5 is the flux passing through the plaquette whose boundary is an oriented parallelogram 
(generalization of the one shown in fig.^ to four dimensions), whose sides are the link piercing 
the 4-cell, and the link along the 1-cell. Such a plaquette has a nonvanishing projection on 
three faces, and has zero projected area on the other three face which are orthogonal to the 
fi direction. Hence, the flux passing through it, is the sum of fluxes passing through the 
faces, which this plaquette has non-vanishing projection, i.e, <^^5 = —Yla^t^^f^a- The two 
commutation relations, Eq. ( p. 11 ) and Eq. (2.12), can be combined into gmg-n = e^^'"" gngm 



by keeping in mind the definition of (75 and <I>^ 
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The reduction of the orbifold lattice action Eq. ( p.5| ) to a single site by using twisted 
boundary condition gives the action of the deformed Type IIB matrix model. The derivation 
of this action is given in the appendix However, the deformed action is easy to understand 
on physical grounds. Here we present a heuristic argument. 

First, note that the lattice action Eq. ( p.5| ) is a sum of a Q-exact and a Q-closed expres- 
sion. Schematically, both types are triangular plaquette terms whose sides are bosonic (B) 
and fermionic (F) Q = 1 supermultiplets. Upon reduction and appropriate field redefinitions, 
all the commutators will be replaced by "deformed commutators", where the deformation 
parameter is the flux passing through the corresponding triangular plaquettes: 

Sciosed ~ Fi[ i3, F2 ] ^ Fi[e^(fl-)i3F2 - e-^^fi^-^FsS] (2.13) 

Since there are p-form supermultiplets which reside on the diagonal of each p-cell, there are 
a number of possible plaquettes. The fluxes associated with each type of plaquette will be 
different. We first write down the deformed type IIB action in manifestly Q = 1 superfields 
and by preserving the U{\)^ global symmetry, and then we explain the plaquette structures. 
The deformed action, which is the main result of this paper, is given by 

Sdef armed = ^ Tr / ( " ^ASgA - A([z^, Z^] + ^ [^'^^^^ Z^,p,]) 

_|_^^g^^^g-«(*PM+*P'^+*CTM+*<T^)/2 -/ii/po- _ gi(*PH+*pf +*<^M+*<T'^)/2 ■^fJ-'^pcr'^ (2 14) 



In this expression, the penultimate and the last line are not separately supersymmetric, but 
their sum add up to a Q-closed form and supersymmetric. The Q = 1 off-shell supersymmetry 
transformations can then be realized in terms of the following superfields, each of which is 
the reduction of the p-form field to a single point 

A = A - eid, 

Z'^ = + v^OiP", z^, 

^IMU £,1^1/ 20i?^jy(z), 

3^5 = Cms - 2eE^5{z) (S'^f^), 
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Z5 = + V2e^'> {Z^upa), ^5 {z^'^n, (2.15) 
where are supersymmetry singlets. The E functions are given by 

F (z\ - p-»*Mi'/2 z z - p+«*M^/2 z z 

E^5{z) = e'^-^^ "^^"/^ z^Z5 - e-'^-^^ ^'"'/^ ^^^^ (2.I6) 

and they give rise to deformed potential of the matrix theory. The E functions are super- 
symmetry singlets since they are functions of z and not Z, and are annihilated by the only 
supersymmetry in the theory: QE = 0. Similarly, the functions E are functions of general 
Q = 1 superfields, and are given by 

E^^ = e"**''^/^ Z'^Z^ - e**''^/^ Z^Z'^ (2.17) 

Notice that, as before, E has an expansion in superspace and the lowest component of E is 
negative hermitian conjug ate of E, i.e; E = -E'^ . 

The rest of this paper is devoted to the study of the deformed action Eq. ( |2.14 ). Even 



though it is not apparent at first glance, it serves as a matrix regularization of the four 
dimensional = 4 SYM theory. This identification requires neither orbifold projection 
nor compactification [22]. We will show how Eq. ( |2.14 ) gives rise to the four dimensional 



= 4 SYM theory, both commutative and non-commutative. In order to do that, it is 
important to understand the structure of the noncommutative moduli space, i.e, identify the 
zero action configuration of the bosonic potential and understand the fluctuations around 
that background. Before moving on we want to give a little bit more detail about the unusual 
plaquette structure and give a visual explanation for the deformations. 
Plaquettes 

The first line of Eq. ( p. 14 ) is the interaction of site superfield A with itself and with 



the link fields. The links traverse 1-cell and the 4-cell diagonals in backward and forward 
directions and vice versa. As it does not surround any area, the commutator is unaltered. 
The second line of Eq. (|2.14 ) is a signed sum over triangular face plaquettes. It is a face- 



link-link (2-cell, 1-cell, 1-cell) plaquette, 2-1-1 for short. The bosonic contribution of this term 
includes the usual twisted Eguchi-Kawai (TEK). This is most easily seen by using the polar 
decomposition of the complex bosonic matrices and by looking to the angular fluctuations. 
As the triangular plaquettes are oppositely oriented, the fluxes comes with opposite signs. 
This is shown in fig.||. 
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po 



Figure 1: Plaquettes types are , respectively 2-1-1, 4-1-3, 3-2-1 and 4-2-2 following the name of the 
p-cell diagonal that the corresponding superfield is residing. The blue link is link fields Z^t, z^, red is 
face diagonal (fermi) multiplet H^j,, green is cube diagonal H"'"^, and black is the diagonal piercing 
the hypercube, Tif^upaiZ^"'"^ 



hypercube-link (or 3-4-1) superfields. The 3-4-1 triangular plaquette has nonvanishing pro- 
jection to only three faces of the hypercube, hence the deformation parameter is just the sum 
of these three fluxes. 

The penultimate line is a triangular plaquette of face-link-cube multiplets ( or 2-1-3) and 
finally the last line is a face- hypercube- face (or 2-4-2). These are shown in the last two figures 
in fig.|l|. The type 2-1-3 plaquette on a fixed Euclidean time-slice is shown in fig.|2[ The 2-1-3 
plaquette in conjunction with 2-1-1 plaquettes provides further insides and is discussed below. 

2.3 An intuitive explanation for the deformations 

Consider a fixed Euclidean time slice on the hypercubic lattice, which is a d = 3 dimensional 
cubic lattice. The background electric flux is parallel to the =fixed hypersurface and 
does not generate a flux for the plaquettes living in that hypersurface. So the deformation 
of the action in X4 = fixed slice will be purely due to magnetic flux <I>jj, where i,j = 1, 2, 3. 
We will show how these fluxes arises by just looking at the orientation and position of the 
plaquettes in figure. 

Let us consider the part of the action that lives on the three dimensional hyperspace. 
The relevant terms in the action are HjjE*-'!^ + H'^*-' [z^, Hjj] where E*-' = [Z*,Z-^]. ^ Both 
are plaquette terms, but of different nature. The first one is the sum of the closed oriented 
loops on the three faces of the cube. The second one involves the diagonal of the cube which 
pierces the volume of the cube, a face diagonal, and an ordinary link. These two types are 
shown in the figure ^. The deformed action is 



^Incidentally, this is also the dimensional reduction of the plaquettes of our rf = 3 lattice action for Q = 8 
target theory. In three dimensional lattice, we make the replacement H*-''' ^^''^X ^-nd X is a supersymmetry 
singlet. See [13]. 



The third line of Eq. ( ^.14 ) is a signed sum of a gauge invariant product of a cube- 
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Figure 2: The constant Euclidean time slice of the hypercubic lattice. All possible plaquettes of type 
2-3-1 and 2-1-1 are made explicit. The flux through a 2-3-1 is the sum of fluxes passing through its 
projection on faces. The color coding is same as -fig.|^. 



The orientation of the gauge invariant product of superfields (shown in figure determine 
the sign of the deformation parameter. The relative signs between 2-3-1 and 2-1-1 type 
plaquettes can also be understood similarly. For example, in figure |2|, any of the 2-3-1 type 
plaquette has nonvanishing projection on two faces and zero projection on the other. However, 
the orientation of the projected 2-3-1 is opposite of the corresponding 2-1-1 plaquettes. This 
accounts for the overall difference of the sign among the fluxes between these two types of 
plaquettes. All the other plaquette terms and deformation parameters in deformed action 
Eq. ( |2.14 ) can be understood similarly. 



3. Noncommutative moduli space 

The classical moduli space of vacua is determined by the vanishing of the bosonic action. The 
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bosonic action of the deformed theory, as expressed in Eq. (|1.1D , is a sum of the d-term and 
E-term. Therefore, the problem at hand reduce to the study of d and £^-flatness, which is a 
set of quadratic matrix polynomial equations ^. 

Surprisingly, the zero action configurations of the deformed Type IIB action turns out 
to be closely related to the TEK action, whose general solutions are well known [7,28]. This 
relation is somewhat unexpected since the the bosonic deformed Type IIB action is expressed 
in terms of ten non-compact, algebra- valued matrices, whereas the TEK is expressed in terms 
of four group valued unitary matrices. The TEK, for a given nonzero deformation parameter, 
does not possess a moduli space, but a discrete set of vacua modulo gauge rotations. This 
discrete set of vacua play an essential role in the description of the moduli space of deformed 
Type IIB. In the rest of this section we first review the zero action configuration of both EK 
and TEK, then explore the moduli space of the deformed Type IIB in connection with TEK. 

3.1 The relation between TEK and deformed Type IIB 

The construction of the deformed Type IIB matrix model, given in appendix follows from 
a dimensional reduction scheme which admits twisted boundary conditions. The reduction by 
using periodic boundary conditions produces the Type IIB matrix model. This is in essence, 
the relation between Eguchi-Kawai (EK) model and twisted Eguchi-Kawai (TEK) model. 

The EK model is a matrix model of d unitary matrices U^, obtained by dimensionally 
reducing the lattice Wilson action to a single point. The action is typically written as 

Sek = -/3 5^ Tr {U^U.Ulul - 1) + h.c. = W^U, - U,U^\'' (3.1) 

Historically, this model had been introduced to work the large Nc limit of the gauge theories 
[29]. It is shown in [29] that the reduced model becomes equivalent to the full theory provided 
the center symmetry (in the reduced model) is not spontaneously broken. However, this does 
not happen to be the case at weak coupling (large /3) [30], and there is indeed a remnant 
of the confinement-deconfinement phase transition at the level of the matrix model. At and 
above the critical value of the coupling constant /3c, the Polyakov loop acquires a vacuum 

®The full examination of the classical moduli space of the arbitrary deformations of Type IIB theory is 
beyond the scope of this paper. The deformations of the matrix model that we constructed are of special type 
and other deformations, which preserve more exact supersymmetry such as Q = 2, 4, 8 and bigger subgroups 
of the global 50(10) R-symmetry, are possible as well. For example, the Fayet-IUiopoulis deformation, mass 
deformation etc. both of which makes the moduli space noncommutative. A classification and analysis of 
the deformations and their consequences can probably be performed along the lines of [9] by applying the 
non-commutative algebraic geometry, and would be interesting. There are also some works exploring the 
deformations of Type IIB model [27] and its bosonic cousins [26] . 
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expectation value. The eigenvalues of the Polyakov loop become nonuniformly distributed, 
and they clump. (For more recent discussions of TEK in context of large Nc equivalences 
see [31], and in the noncommutative gauge theory context see [32]) 

The twisted EK is a variant of the EK model, which is introduced to prevent spontaneous 
symmetry breaking, hence is a cure for the clumping of eigenvalues. (Another solution is 
quenching [30] ) The action is altered in such a way that the eigenvalues of the unitary matrices 
are uniformly distributed even at the very weak coupling. The TEK action is given by [2] 

Stek = -/5^TT{e-'''^^U^U,Ulul - 1) + h.c. = J^Tr le-'^-^/^U^.U, - e^^-Z^^.C/^l^ (3.2) 



where 6**^" is the flux factor as in Eq. pTll ). Notice that the TEK action can be thought 



part of the full deformed Type IIB action. This can be done by doing a polar decomposition of 
the complex matrices as = H/^iU/^i and looking to a configuration for which the radial modes 
are proportional to identity, as in [33]. Freezing the radial fluctuations and paying attention 
to only angular ones yields Ylfij^u ~ Stek-, where £"s are given in Eq. ( |2.17|) . 

The twisting (deformation) in Eq. ( p.2| ) is an element of Zjv^ and has an important 
implication. It forces the eigenvalues of the Polyakov loops to be uniformly distributed even 
at weak coupling, preventing the Polyakov loops from acquiring vacuum expectation values. 
The zero action configurations of TEK action are given by 

C70C/0 = e^*M.f70f70 (3 3) 

which has well-known twist-eating solutions. These solutions are expressed in terms of clock 
Pl and shift Ql matrices given by {PL)jk = e^'^^'^/^5jk and {QL)jk = Sj+i^k- 

We consider the simplest deformations in U{Nc) matrix theory corresponding to n units 
of electric and magnetic flux in z-direction. The flux matrix comes into the canonical form: 

27rn 

^t^u = —l2^i(J2 (3.4) 

where I2 is two dimensional identity and (T2 is the antisymmetric Pauli matrix. The zero 
action configuration for Nc = and n = L background can be written as 

U^ = Pl^1l, U^ = Ql®Il. UI = 1l^Pl, UI = Il®Ql (3.5) 

The solution is unique modulo gauge transformations. Its significance for the deformed type 
IIB will be discussed in depth in the next section, but at this point it is useful to realize that 
these are generators of a four dimensional (fuzzy) torus. 
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3.2 The noncommutative moduli space of deformed Type IIB 

The zero action configurations of the Eq. ( 2.14| ) are determined by the d and E fiatness 
conditions 

5 

id= = 0, 

m=l 

^^5 = z-z^ - e-^^'-^M z^z^ = (3.6) 

which are a total of eleven quadratic matrix polynomial equations. The d-flatness condition is 
unaltered with respect to undeformed case. The matrices satisfying the S-flatness condition 
will also satisfy the d-term condition, following the analysis of [34] . Henceforth, our goal is 
to obtain the solution of the E = {). 

Let us first recall the moduli space of the (undeformed) Type IIB theory, given by d (as 
above) and iiJ-flatness conditions = [z™,^"] = 0. Solutions are mutually commuting, 
diagonal matrices with complex eigenvalues. Hence, the moduli space of the matrix model 
Eq. ( ^^ ) with gauge group U{Nc) is 

M = C^^'/Sn, (3.7) 

where Sn^ is the Weyl group of U{Nc). This moduli space is a commutative one since the 
matrices are simultaneously diagonalizable. 

We now describe the noncommutative moduli space, given by Eq. (|3.6| ). The equation 
in the penultimate line of Eq. ( p.6| ) is, in fact, the same as the equation for the minima of 
the twisted EK model Eq. ( |3.3| ), with one essential difference. In TEK, the matrices are 
group valued, unitary matrices. Whereas, the z^ matrices entering Eq. ( ^.6| ) are noncompact 
algebra- valued, complex matrices. The solutions of Eq. ( p.3| ) will satisfy the Eq. ( p.6| ) as well. 
Then, the solution of Eq. ( p. 61 ) will be of the form 

Zfi = Cfj^U^, no sum (3-8) 

where is a complex moduli which is associated with the rigid rotations and overall scaling 
of the eigenvalues. The equation = will be satisfied if 

Z5 = c^U^ = C5 (^llU^^ , (3.9) 

where C5 is an independent complex number. The 25, apart from the scale factor C5, is not 
an independent matrix and its form is fixed by the other four C/^'s. 
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a) 

Figure 3: The distribution of eigenvalues of one of the matrices. The moduli fields corresponds 
to the rigid rotations of eigenvalues and their overall scaling, a) The k — 1 case. The target theory 
has a U{1) gauge group, b) The k — 2 case. For coinciding c^^ and c^, there is an enhanced gauge 
symmetry group U{2), otherwise the gauge group is U{1) x U{1). 

The space of the minima of the deformed type IIB action is therefore parametrized by five 
continuous parameters, ci, . . . , C5 and a finite set of unitary matrices satisfying the constraint 
Eq. ( |3.3D . The eigenvalues of the twist-eating matrices are roots of unity. The eigenvalues 
of each Zm are uniformly distributed over a circle of radius c^, see fig.|3| a). The modulus and 
the phase of Cm correspond to the overall scaling and the rigid rotation of the eigenvalues, 
as noted earlier. To be more precise, for Nc = L"^, and flux n = L, the eigenvalues of the 
matrix Zm are partitioned into L groups, located at Cme^^^^^^ where A; = 0, ... L — 1. For our 
interpretation of the deformed matrix model as a four dimensional field theory, these moduli 
modes are essential ^. 

Modulo these five complex numbers, and an extra unitary matrix (which is fixed by 
the other four U^), the set of zero action configurations of the deformed Type IIB matrix 
model Eq. ( |3.6D is same as the one of the corresponding TEK (with the same background 
flux) . We can express the moduli space of the deformed theory as 

Mdeformed = X MtEK (3.10) 

Here Mtek corresponds to a finite, discrete set of vacua which happens to coincide with the 



minima of corresponding TEK Eq. (p.sp As noted in section 3.1, the procedure to construct 



Mtek is discussed in depth in literature and books (see for example [35]). 



^The moduli space shown in figure ^ was raised by Nima Arkani-Hamed almost a year ago in a discussion. 
I thank him for sharing his ideas with me. 
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The deformation changes the structure of the classical moduli space drastically ^ , turning 
the classical commutative moduli space A4 = C^^^/Snc in which the eigenvalues can take 
any value (without increasing the action) to a space in which there are only five free complex 
parameters. Recalling the polar decomposition of a complex GL(Nc,C) into a product of 
a Hermitian positive definite matrix and a unitary matrix, = H^U^, we observe that, 
in the moduli space of the deformed theory, the radial coordinates are forced to be equal: 

= Ic^IIat^. The nonuniformities disrespecting the Z^, in the eigenvalue distribution, in 
either radial or angular direction, increase the action. 

The reader may also wonder if there are other solutions to Eq. (|3.6| ) and if any, what 
do they correspond to. There are indeed such solutions and the ones for which at least two 
out of five Zm are zero are of different character. For example, in = 2^ = case, it is not 



possible to generate a four dimensional spacetime. As we will see in section iA, the main 
moral of Eq. ( |3.8| ) is that it forms a complete basis for GL{Nc,C) matrices, which in turn 
becomes the basis for the emergent lattice theory. The case for which at least two out of five 
matrices is zero a full basis for GL{Nc, C) can not be constructed and we do not obtain a four 
dimensional target theory. However, it is still possible to obtain a two dimensional target 
theory. If at least four out of five matrices are zero, then it is not even possible to obtain a 
two dimensional theory. The technique for the study of the full set of solutions is given in 
ref . [9] , and a similar deformation of the superpotential of = 4 is worked out there as an 
example. We refer the reader to ref. [9] for further discussions. 
(Ir)reducihle Representations 

Up to this point, we have discussed the theory with = and n = L, where ^>^,y = 
(E> i<72- As we will see momentarily, this background supports a U{1) target theory on a 
four dimensional lattice with sites. In order to get a U (k) theory on the same lattice, all 
we need to do is to take Nc = L'^k and n = Lk. In this case, the moduli space is described 
by matrices 

z"' = Cm'^U^, where = Diag (c^, . . . , c^) (3.11) 

where Cm is diagonal k x k moduli matrix and the solution can be thought as a direct sum 
of irreducible solutions in the form 0*L]^ c^f^m = Cm U^. Therefore, the moduli space can 
be expressed as 

Mde for med = C^V^fc ® MtEK (3.12) 

*For special values of flux, such as ^^i, = except for <l>i2, the moduli space becomes (locally) C^^"" x 
X Mtek2 i.e, the eigenvalues in the 2:3 and Z4, start to move freely again and Mtek2 is the zero action 
configuration of two dimensional TEK. This deformation can be shown to be relevant to construct the TV = 
(8, 8) theory in two dimensions. Finally, when all the deformations are turned off, we recover M = C^^" /Sn^ 
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The four dimensional lattice basis in terms of zero action configurations is encoded into 
Mtek, and will be given explicitly in the next section. We will see how the noncommutativity 
of the moduli space is tied to the noncommutativity of lattice base space. 

Before moving on, we want to make one more observation. Geometrically, the k = 1 
describes the eigenvalues residing on the surface of a four dimensional (fuzzy) hypertorus T^. 
For k > 2, there are multiple, concentric fuzzy tori T^, see figj^). When these concentric 
tori are on top of each other, there is an enhanced gauge symmetry group U{k). If all of them 
are separated, the gauge group of the lattice theory is U{1)^. This phenomena is same as the 
emergence of U {k) gauge theory on a stack of k D-branes wrapped on a hypertorus. We will 
comment more on that in section ^. 



The emergence of (noncommutative) base space 



4.1 Fluctuations 



In this section, we show that the deformed Type IIB action Eq. ( 2.14| ), without any orbifold 



constraints whatsoever, is a lattice regularization of = 4 SYM theory. The emergent lattice 
is noncommutative. However, depending on the deformation parameter, it is possible to find 
two continuum limits: one is regular = 4 SYM, and the other is its noncommutative 
counterpart. 

Let us start with the deformed Type IIB action Eq. (|2.14 ) with gauge group U{Nc) 
where Nc = L'^k and with the flux matrix (^^^ = ^(l2 ® '1(^2)- In section 3.2, we found the 



zero action configuration of deformed type IIB. We will analyze the bosonic and fermionic 
fluctuations around this background. To do that, we decompose the bosonic matrices into 
fluctuations and background as 

= ZmU^^ = no sum (4.1) 

where is a background field configuration and are fluctuations around it. The change 
in naming of the background is due to the physical interpretation, that the will have 
momentarily, as displacement (discrete translation) operator on the lattice. Substituting 
Eq. ( [4.1] ) into the bosonic part of the deformed action Eq. ( |2.14 ) yields 

5 



Sb = ^ 



m=l m,n 

where the phase factors disappeared upon the use of the Eq. ( [3.31) . 



(4.2) 
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The fluctuations entering the action are U{Nc) algebra valued and are conventionally 
expressed as z = z ^t^; A = l,...N^; where are Hermitian generators and satisfy the 
standard commutation relations [t'^jf'^] = if^^Qt'~" . Instead, we will use a basis which 
makes the mapping of the Nc x Nc (where Nc = Li^k) matrix into U{k) fields on lattice 
clearer. We introduce Jp = e-^'^(PiP2+P3P4)/L(f/0)pi(f/0)P2(f/0)P3([/0)P4^ where p^, = 1,...L 
which form a complete and orthogonal set for x matrices, GL(L^,C). The Jp basis is 
identifiable as the dual momentum basis for the lattice and p is an integer vector labeling 
points in Brillouin zone, By a Fourier transform, we construct the coordinate basis An, 
where 

^" = 14 E '^pe'"*^-"/'' (4.3) 

Hence, we can express the map between a z fluctuation matrix and fields in the dual 
momentum space lattice (or equivalently real space lattice ) as 

z = Y, ^^t^ = E ® 4k =Y.^p®4 I= ^ ® An (4.4) 

We are essentially expressing a GL{L'^k,C) matrix by using complementary complete basis 
matrices, Jp G GL(L^, C) and G GL{k, C). Notice that this representation makes clear how 
color degrees of freedom of the deformed matrix theory turns into a noncommutative space 
time in the lattice theory. The properties of the basis matrices are given in [35]. The beauty 
of this basis is that the zero-action configurations = acts as displacement operator on 
the position space lattice. Namely, 

Df, (An Zn) DI = An-e^ ^n (4.5) 

This relation provides a nearest neighbor interaction on the lattice. By using the orthogonality 
of basis matrices Jp and Eq. ( |4.5|) , the bosonic action can be expressed as 



nezt 



2^ 

Cy ^ (^Tt.n— Bm * ^m,n~em ^m,n * ^m,n)) + ^ ^ | ( 



_ |2 
Zm. ri ^ Zn rt4.p — 171 < > Tl\ 



m,n 



(4.6) 



This is indeed just the noncommutative generalization of bosonic part of Eq. (^]^). Performing 
the same analysis at the level of superfields, it is easy to show that the deformed Type IIB 



action Eq. ( 2.14 ) may be rewritten as 
Tr , f / 1 

S = J dOi— 2 * dgAn — An * (^m,n-em * '^m,n~e^ — Z^.n * ^r; 
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The lattice star product is defined by 

which is a non-local product and 9' = 2/L is a dimensionless noncommutativity parameter 
on the lattice. 

4.2 Classical spectrum 

The action Eq. (4.7) has no explicit lattice spacing, but it has a classical moduli space for 



which the bosonic lattice action vanishes. The distance from the origin of the moduli space is 
identified as inverse lattice spacing, and hence the continuum limit is defined as a trajectory 
out to infinity in the moduli space. 

The zero action configuration of the lattice corresponds to site independent diagonal 
matrices Zm,n = Cm = Diag (c^, . . . , c^), hence the moduli space is given by = C^^ /Sk- 
Notice that this is (and has to be) the moduli space of the deformed Type IIB action, as 
Eq. ( [4.7] ) is merely a rewriting. 

We expand the action around a U{k) gauge symmetry preserving, hypercubic spacetime 
lattice point: 

zi^ = Zf,,n = 4 = ^5,n = (4.9) 

\/2a 



We use a cartesian decomposition of the bosonic matrices around the point Eq. ( [4.9[ ) as 



where S and V are identified as six scalars and four gauge bosons of the target = 4 
supersymmetric Yang-Mills theory. 

The classical spectrum can be found by diagonalizing the kinetic terms, similar to the 
appendix B in [4]. This can be done by introducing the Fourier transform of a field \I'n 

peB.z. 
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where the = "^p^ G [— 7r/a,7r/a) is discrete momenta in hypercubic Brillouin zone (B.Z.), 
and Pfj_ is integer in the range [— L/2, L/2). This diagonalize the kinetic terms for the bosonic 
action 

^ E [E^'^.P^P^«.-P + ^M,pG^^,pK,-p] (4.12) 

^ peB.Z. a=l 

where we define 

9 4 
EE - sin ^ , Ml = Y,^l, G^,,p = M^J^, - V^V,e-"'^^P--P'^^'^ (4.13) 



The six scalars are degenerate for a given momenta and have eigenvalue Mp. The eigenvalues 
for the gauge boson matrix G^i/,p are a three fold degenerate Mp and zero. The zero mode of 
G is the consequence of gauge invariance. The essential point here is the absence of doublers 
in the bosonic spectrum, Mp, which is only zero at the origin of the Brillouin zone and which 
reaches its maximum at the corners. The exact supersymmetry implies there are no fermion 
doublers either (see [12] for an explicit analysis.) The absence of the fermionic doublers is 
also a consequence of associating the sixteen fermions with the p-cells of the lattice (see [24]). 
^ . In the continuum limit, which is defined by 

a ^ 0, L ^ oo, La = fixed, (4.14) 

we obtain — > and Mp p^ and G/ij/,p — > ip^^^iu — Pfj.Pi/) correctly reproducing the 
kinetic terms of bosonic action. 
Moduli space: A second pass 

We reconsider the eigenvalue distribution shown in fig.^ and make various observations. 
For simplicity, consider the case k = 1. In finding the spectrum, we expand the fields around 
the point Eq. ( [4.S| ) in moduli space. This means setting = ^ and C5 = 0, where a is 
interpreted as the lattice spacing. Turning on C5 is related to deforming the hypercubic lattice 
to the lattice and tilting the T^, and we ignore this possibility for simplicity. Consider the 
eigenvalue distribution on the circle shown in fig.^, which is one of the directions on fuzzy 
on which eigenvalues are distributed. The circumference of the circle is 

27r 

27rc^ = — (4.15) 
a 



^Notice that the spectrum Eq. (4.13) is not the same as the spectrum of A4 lattice, given in [4]. This is 
due to difference in the lattice structures, ( Al versus hypercubic ), and consequently their Brillouin zones 
{A4 versus hypercubic). This is analogous to the difference in the acoustic phonon spectrum of body-centered 
cubic and simple cubic lattices in solid state physics. However, in our case, the long-wavelength hmit of these 
two lattice theories coincides and become the — 4 in four dimensions. 
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which is, in fact, the extent of the Brihouin zone in the ^ direction. The separation between 
two nearest-neighbor in moduh space along the circle (geodesic distance on the torus) is 
given by 2itc^/L = |^ and this corresponds to the minimal lattice momenta. The geodesic 
separation between two eigenvalues separated by S [— L/2,L/2) units is given by = 
"WaPi^ S [— 7r/a,7r/a) and is the discrete momenta in Brillouin zone of hypercubic lattice. As 
a last remark, notice that the external distance between the same two points as above is 
|sin^y^, the acoustic phonon spectrum in one dimension. The continuum limit is taken as 
a trajectory out to infinity in moduli space, by keeping the volume and four dimensional 
coupling constant fixed. Hence, the circumference of eigenvalue circle in moduli space (and 
extent of the Brillouin zone) becomes G [— vr/a, vr/a) [— cxd,cxd). It seems natural to 
interpret the four dimensional torus that the eigenvalues reside as the Brillouin zone, and the 
action as the momentum space action, with momenta summed over the Brillouin zone. 

In the derivation of Eq. ( [4.7] ) , we used a Fourier transform, which took us from momentum 
space representation (Jp) of the lattice to the position space representation (An). Hence, the 
inverse extend of Brillouin zone (divided by 2ti ) becomes the UV cutoff. The inverse of 
minimal momenta (divided by 2it) becomes the IR cutoff {La), the size of the spacetime 
torus. The four dimensional target field theory lives on this spacetime torus, whose extent 
remains fixed as we take the continuum limit. It should be emphasized that the target space 
is not the space that the eigenvalues reside, but its Fourier transform. 

4.3 Localization of interactions 

In this subsection, we will describe the circumstances under which the non-locality in Eq. ( ^.Tj ) 
disappears and the interaction becomes local in spacetime. Upon that, the continuum limit 



of the action becomes the same as the continuum limit of Eq. (2.5) and reproduces the N = A 
SYM in four dimensions. Since the calculations are very similar to [4], they will not be 
duplicated here. At the end, we also comment briefly on the non-commutative spacetime 
limit. 



The interaction terms of the lattice action Eq. (4.7) are non-local, involving a sum over 



all lattice points. In the continuum limit, the star-product becomes 

-^i{\v)i^-^2{^) = J d^x d'^y ^i{x)K{x-w,y -w)^2iy) ■ (4.16) 
The kernel of the integral is 

K(yL - w, y - W) = g-2i(x-w)^(0-i)M-(y-w), 

vr^ det @ 
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and the dimensionful noncommutativity parameter is given by 



2-K 



e^,^ = — (I2 X ia2)^u- (4.18) 



Let us consider the continuum hmits in which the volume of the target space is kept 
fixed. The length scale associated with noncommutativity is the square root of Q12, denoted 
as \/0 for short. Then, for y/Q/{La) 0, i.e, the limit at which noncommutativity goes to 
zero, the kernel localizes sharply and becomes a Dirac delta function: 

K(x - w,y - w) ^ (5(x - w)J(y - w), (4.19) 

yielding a local product out of the star product. From the spacetime point of view, trading 
the nonlocality of the product of functions back to spacetime noncommutativity, this limit 
yields a commutative spacetime. This is certainly an interesting limit as it reproduces the 
ordinary = 4 SYM target theory. 

Explicit calculation shows that, in the continuum limit, we obtain the action of = 4 
SYM at tree level. 

S = ^ y d^xTr [^V^uVf,, + ^{Df^Saf - \ [Sa,Sbf + C {D^t^u; + iV^+a [Sa, w])) (4.20) 

where = + i[V^, ■ ] is the covariant derivative and V^u = —i[Dj_i, D^] is the nonabelian 
field strength. The 5*0(10) gamma matrices are basically a reshuffling of the gamma 
matrix basis introduced in Eq. (^). They are given as = T^, T4+m = Eio = T^, 

following the conventions of section 2.2 of [4]. The charge conjugation matrix C is unchanged 



with respect to Eq. (2.1 



Even though the deformed type IIB theory Eq. ( p.l4| ) , hence lattice theory Eq. ( [4.7| ), 
possesses only Q = 1 exact supersymmetry and a U{1)^ global symmetry, the target theory 
has full Q = 16 supersymmetries and manifest Lorentz invariance along with a global S0{6) 
R-symmetry. 

We believe it is also possible to obtain a lattice action and its continuum for the product 
gauge group YliU{ki) where J2if^i = k hy expanding around a point in the moduli space 
given by (zn) = "^^Ifc + (©jCjlfci)- In the continuum limit, we take the lattice spacing to 
zero by keeping q's fixed, which produce a YliU{ki) gauge theory. It is desirable to know 
whether one can obtain a lattice or matrix regularization for orthogonal and symplectic gauge 
groups by either performing an orientifold projection on the orbifold lattice or (deformed) 
type IIB matrix model (The symplectic type IIB is discussed in [36]). Such constructions 
may be useful in nonperturbative study of electric-magnetic duality [37]. 
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Noncommutative continuum limit 

We believe it is also possible to obtain the noncommutative continuum = 4 theory 
with the appropriate choices of the background flux matrix ^^u, and hence appropriate zero 
action configurations U^, for example the ones given in [16]. In order to achieve this limit, 
one needs to take 0' fixed (hence fixed) along with a fixed finite volume. In this case, there 
is a nonlocality of order \/0 in the star product of the continuum theory. At distances longer 
then V&, the interactions effectively shuts off because of the fast oscillation of the kernel. 
The action of the noncommutative = 4 SYM target theory can be obtained by replacing 
the product of fields in Eq. ( [4.20 ) by nonlocal continuum star product given in Eq. ( 4.16 ). 



5. Collective excitations of D(-l) branes as D3 branes 

The results obtained within this paper can be interpreted within the string theory framework. 
In order to do that, we need to specify the world- volume theories of the corresponding D- 
brane configurations. The type IIB matrix theory Eq. ( p.l| ) describes the dynamics of a 
cluster of D(-l)-branes (which are sometimes referred as D-instantons). These are points in 
ten dimensional Euclidean space. In this sense a Dp-brane is a p+1 (Euclidean) dimensional 
hypersurface. The main result of this paper may be interpreted as follows: Upon deforming 
the type IIB matrix theory, a particular collective excitation of D(-l)-branes becomes the 
classical vacuum state of the theory. Essentially, the D(-l) branes, instead of being at some 
random points on the moduli space of the undeformed theory, lie on the surface of a four 
dimensional torus T^. Analyzing the fiuctuations around this D(-l) background yields AA = 4 
supersymmetric Yang-Mills theory, which is the world-volume theory of a collection of D3- 
branes. Hence the D3-brane, which is not present in the fundamental description of the 
theory, may be considered as an emergent brane. In the rest of this section, we will describe 
these statement in a little bit more detail. 

The type IIB action Eq. ( p.lj ) (or Eq. (^)) with gauge group U{Nc) describes a collection 
of Nc D(-l)-branes (see for example [38]). This theory possess Q = 16 supersymmetries and 
a global 50(10) symmetry. In the latter form of the action Eq. (|2.8| ), only the SU{5) x 
C/(l) subgroup of 50(10) is manifest. The commutative moduli space is (locally) R^'^^'^ (or 
equivalently C^^'^). The eigenvalues of these matrices {va or z"^) are identified as position of 
D(-l) branes. In this framework, the D-branes are pointlike objects on the moduli space. 

The deformed Type IIB action Eq. ( |2.14| ) is a supersymmetric deformation respectful to 
a global U{1)^ symmetry and Q = 1 supersymmetry. The deformation renders the moduli 
space (defined by E and d-fiatness conditions given in Eq. ( ^.6]) ) a noncommutative one, and 
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forces the eigenvalues of Zm matrices (the D(-l) branes) to be uniformly distributed on a circle 
as in^. However, as the matrices zi, . . . are not simultaneously diagonalizable, the positions 
of the D(-l) branes are indeterminate and they lie on the surface of a fuzzy four dimensional 
torus. The classical zero action configuration, given in Eq. ( |3.8| ) and Eq. (|3.9| ), is a collective 
excitation of D(-l)-branes and is responsible for the emergence of higher dimensional branes 
out of the deformed Type IIB matrix theory. Also notice that, the base space of the target 
theory is not the space that the eigenvalues reside, but the dual space obtained by Fourier 



transform Eq. (^j) in section iA. The continuum limit is taken by keeping the size of this 
dual space fixed. This is similar to the Seiberg-Witten limit where the open string metric and 
coupling constant rather then the closed string parameters are kept fixed [39] (also see [10,40]). 

It is also reasonable to ask which supergravity field or type IIB string theory field may 
be responsible for this excitation of D(-l) branes. To answer this question, let us first recall 
that the world-volume of a Dp-brane couples to a (p+l) form potential A^'^^^\ given in the 
form A(P+i),or explicitly as / ^'^i-'^p+i ^ . . . J^e^i-*f+i(iP+V, where are world- 
sheet coordinates and the integration is over the world-sheet of Dp-brane. The D(-l) branes, 
whose world-sheet is zero dimensional, couples to a zero form. The configuration we consider 
above carries N^. units of D(-l)-brane charge. Moreover, the fact that the matrices z^ do not 
commute with each other can be interpreted as the existence of higher dimensional branes in 
the system (see [22] for a review). In our case, the emergent higher dimensional branes are k 
D3-branes wrapped on a four dimensional spacetime hypertorus, where k is associated with 
the rank of the gauge group U {k) of the target theory. It is, therefore, tempting to think that 
the deformations introduced in Eq. ( p. 14 ) on the matrix theory side corresponds to turning 



on the RR five-form self dual field strength F^^^ on the type IIB string theory. The F^^^ arises 
from a four-form field A^^\ which, as noted above, couples to D3-branes. The emergence of 
D3 branes can be thought as a special case of the dielectric brane eff'ect introduced by Myers, 
which is analogous to polarization of a neutral atoms in an electric field [41]. Similarly, the 
collective state of D(-l) branes can be thought of as analogous to the emergence of a dipole 
upon the application of F^^^ field. However, the argument we have given here is incomplete as 
it does not specify the geometry or topology of the emergent object. For example, if we were 
to construct the analog oi M = 1* theory in four dimension [42], i.e, a mass deformed matrix 
theory [43], the emergent brane would be a (euclidean) Dl brane wrapped on a two sphere 
5^, whereas the deformation discussed in footnote 8 gives a Dl wrapped on a two torus T^. 
Both of these can be thought as polarization of D(-l) branes. The difference is similar to 
the emergence of dipole or quadrupole moments upon the application of appropriate electric 
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fields. Certainly, the string theory realization of the deformations warrants further study. 

It is also possible to recover the known membranes constructions appearing in Type IIB 
matrix theory from the large Nc limits of our formulation. Let us recall the zero action 
configurations = c^U^^^ and restrict our attention to the small fluctuations of the 
and small ^^u- Then the conditions describing the vacua uj^^^ui^^ = e^^'^'^ui^^U^^ become 
= —i^^y, which describes (locally) a planar membrane, where are the fluctu- 
ations. The latter equation does not admit a solution in terms of finite matrices, and is 
usually considered to be associated with the existence of a D(p-|-2) in a system of Dp branes. 
However, notice that such a configuration is not stable within the Type IIB, but its decay is 
is suppressed in the large [23]. The zero action configurations of the deformed type IIB, 
on the other hand, admit solutions in terms of finite dimensional matrices, and, more impor- 
tantly, are stable (This point is also emphasized by Dorey [10]). It is also worth mentioning 
that the analogs of a planar membrane, in the notation of Eq. ( p. 14 ), is a Fayet-Illiopoulis 



term, given in the form of a d-term deformation CTr / d9 A = —(d, where C is the deforma- 
tion parameter. This alters the d-flatness condition into Ylm=i['^rn, z"^] — C = which only 
admits solution in terms of infinite dimensional matrices. 
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A. Derivation of deformed type IIB action 



In this appendix we show that the dimensional reduction of the lattice action Eq. (2.5) to 



one site model, by using the twisted boundary conditions Eq. (2.10), gives the deformed type 



IIB action Eq. ( 2.14 ). We will first perform the reduction. However, this will produce some 
extraneous gauge rotation matrices, which we will get rid of by appropriate field redefinitions. 
The first two term in action Eq. (2.5) reduce to a form 



Tr 



-^AdgA - A[{gl^Zm), (Z^ffm)] 



(A.l) 



which has no flux terms. This is due to the fact that the interaction of site fields with the 
link fields does not surround an area. 
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Next, we consider the interactions of fermi multiplets H^n with holomorphic E™" func- 
tion. First, notice that the 9 component of the Smn fermi multiplet is a composite field, and 
upon reduction it becomes 

= gn9m{glnZm){9lzn) - 9m9n{9liZn){9L'^m) = 9m9n (e"**™" (^m^m) (^n^n) " (fi^n^n) (ffL2m))(A.2) 

where we used 't Hooft algebra. Hence the third term in the action becomes 



Tr xprnni rp 



{9i9lSmn)[ (Z™9^)(Z"5n) " 6**™" (Z^ffJ (Z^^n^)] „ (A.3) 



Even though we combined the two terms H^^^E'^^ + 2H^5E'^^ into SmnE™'"', these are plaque- 
ttes of different types, namely a 2-1-1 and 4-1-3 types (see fig.|l]). 

We are done with the Q-exact part of the lattice action. The Q-closed part of action 
includes plaquettes of two types. One is 4-2-2 type, shown in fig.|| and the other is 3-1-2 
type shown in fig.^ and more explicitly, in fig.^. The first type involves triangular loops 
composed of the 4-cell link field and two 2-cell fermionic links. The second type involves a 
3-cell fermionic link, a 2-cell fermionic link and a bosonic 1-cell link field. Let us first consider 
the first type and its reduction. Using the action Eq. ( |2.5D and setting p = 5, we obtain 



\/2 

— ^ d^[y^n(^5,n— e5"p(T,n+ef,.+ei, ~ '^po-,n— Gp— Go- ^5,n+e^+e^ ) ^ 

a/2 

—e^"''"'Sf,^{{glz5g5){gf,g^Sp„glgl) - {glgl'Bp„gpg„){g^g^Z5glgl)) 

_gM^P<7(^^tgtH^^) ((55^5) {959fi9u9p9cT) (gUl^pa) - {9l9l'^pa)(9p9a9ii9u9^){9\M) (^-4) 

In this expression, since the combination {g^gp9u9p9(T) = Rfiupa appears with e'^'^f"^ tensor, 
with the use of the definition of g^ Eq. ( p. 10 ), we see that it is proportional to the identity 



matrix and can be pulled out of the parenthesis. By using the 't Hooft's algebra Eq. (2.11) 



in the reordering of the second term, we express Eq. (A. 4) in a simpler form as 



"^e'^'^P^Rp.pMgl^pu) {{gh) {gWpSp^) - e^(*'"'+*-+*-+*-^)(5i5^H,.)(4s5)j (A.5) 
where the exponent is a sum of the fluxes passing through all surfaces but fiv and pa. 



The 3-1-2 type plaquettes in Sdosed can be obtained by setting p / 5 in Eq. (2.5) 
Following very similar steps to the reduction performed above, we obtain 
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V2 



2 e^^'^^'^R.^pMoi'^,,^) ({gpp) (5ki3.5)e 



llgl^^Mp-^p)) (A.6) 



By doing a field redefinition, we can easily remove all the gauge rotation matrices. An 
appropriate choice that we can read off from above considerations is: 



Zr. 



ynym^rnn 



(A.7) 



The action of the deformed theory takes the simple form: 
Tr 

■_^iJ.upan c 



_ 1 

d 9 

9^ 



S, 
V2 



^AdeA - A[zm, Z^] + ^H^„(Z-Z- - e^*-Z"Z™) 



+ 



'^p.vpa^pv 



per ^ ^pa^bj 



(A.8) 



where the usual commutators are altered to commutators involving fluxes. 

However, the above form of the action is not completely satisfactory. First, there are 
phases Rpupa in second line of Eq. ( A.8 ), and second, the redefined variables are not perfectly 
antisymmetric. For example, the redefined H^j^ 7^ —'^up, but S^,^ = — e~**'"'Hi,^. In fact, 
our redefinition of the fields does not respect to the antisymmetry of the p-form fields (or it 
only respects up to a phase factor). These two problems are indeed tied, and can be cured 
by a better field redefinition which makes the antisymmetry properties manifest: 

yo-ypyu^ 



■'vpcj 



K^ap+^au+^ap+^pu+^pp+^^p)/^ n n n 7 

yo-ypyi/ypz 



-pupa _^ jpiypa- 



{AS 



By using the fact that the combination R 



pupa ' 



-i{^p.+^^p+^^a+^.p+^u.+^p.)/2 jg ^ constant 



phase, and some simple algebra, we obtain the action Eq. (2.14). This action, as already 
stated, is a Q = 1 supersymmetry, and U^Vf global symmetry preserving deformation of the 
type IIB matrix theory. 
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